Large bulk resistivity and surface quantum oscillations in the topological insulator 
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Topological insulators are predicted to present novel surface transport phenomena, but their 
experimental studies have been hindered by a metallic bulk conduction that overwhelms the surface 
transport. We show that a new topological insulator, Bi2Te2Se, presents a high resistivity exceeding 
1 Qcm and a variable-range hopping behavior, and yet presents Shubnikov-de Haas oscillations 
coming from the surface Dirac fermions. Furthermore, we have been able to clarify both the bulk 
and surface transport channels, establishing a comprehensive understanding of the transport in this 
material. Our results demonstrate that Bi2Te2Se is the best material to date for studying the surface 
quantum transport in a topological insulator. 

PACS numbers: 73.25.+i, 71.18.+y, 73.20.At, 72.20.My, 71.55.Ht 



The three-dimensional (3D) topological insulator (TI) 
is characterized by a non-trivial topology^ of the 
bulk wave function, and it represents a new topological 
quantum state realized in a band insulator. In theory, 
3D TIs are insulating in the bulk and unusual metallic 
surface states consisting of spin-filtered Dirac fermions 
give rise to interesting surface transport phenomenal— 
In reality, however, TI samples available today are in- 
variably conducting in the bulk, and charge transport is 
always dominated by the bulk current £r— Therefore, to 
exploit the novel surface transport properties of topolog- 
ical insulators, it is crucial to achieve a bulk-insulating 
state in a TI material. 

Among the recently discovered TIs, Bi2Se3 has been 
the most attractive because of its simple surface-state 
structure. 3 Unfortunately, near-stoichiometric Bi2Se3 is 
always n-type owing to a large amount of Se vacan- 
cies. An isostructural material Bi2Te3 can be grown 
as p-type^S but usually it is also highly metallic, most 
likely due to anti-site defects which are promoted by 
close electronegativities of Bi and Te. Significant efforts 
have been made£~— to achieve bulk insulating behav- 
ior in Bi2Se3 and Bi2Te3; however, while an increase 
in resistivity with decreasing temperature has been ob- 
served, so far the bulk remains to be essentially a metal 
and a clearly insulating temperature dependence, such as 
the variable-range hopping (VRH) behavior,— has never 
been reported. For example, by growing a Bi2Te3 sin- 
gle crystal with a compositional gradient, it was possible 
to observe a resistivity upturn at low temperature and 
to measure the surface quantum transport,— but the re- 
sistivity remained low (< 12 mficm) in absolute terms 
and the surface contribution to the transport did not ex- 
ceed ~0.3%»i£ The situation is essentially the sained— in 
Bi 2 Se 3 . 

Given this difficulty, searching for a new TI material 
better suited for achieving the bulk insulating state is 
obviously important. In this paper, we report that a new 
TI material, Bi2Te2Se, which has an ordered tetradymite 
structure 15 [Fig. 1(a)] with the basic quintuple-layer unit 
of Te-Bi-Se-Bi-Te and was recently confirmed to have a 



topological surface stated has desirable characteristics 
for surface transport studies. We found that high-quality 
single crystals of ordered Bi2Te2Se show a high resistivity 
exceeding 1 ficm, together with a variable-range hopping 
(VRH) behavior which is a hallmark of an insulator; yet, 
it presents Shubnikov-de Haas (SdH) oscillations which 
signify the 2D surface state consistent with the topo- 
logical one observed by the angle-resolved photoemission 
spectroscopy (ARPES)^ 

By examining the difference in the doping chemistry 
between Bi2Se3 and Bi2Te3, one may understand that 
the ordered Bi2Te2Se has reasons to be superior: i) The 
formation of Se vacancies is expected to be suppressed, 
because the Se trapped between two Bi atoms is less 
exposed to evaporation due to stronger chemical bond- 
ing with Bi in this position; ii) The formation of the 
anti-site defects between Te and Bi is also expected to 
be suppressed because of preferable bonding between Se 
and Bi in contrast to Se-Te bonding; iii) Ordered nature 
minimizes the additional disorder that could be caused 
by Se/Te randomness. In this work, single crystals of 
Bi2Te2Se were grown by melting high purity (6N) ele- 
ments of Bi, Te and Se with a molar ratio of 2:1.95:1.05 
at 850°C for two days in evacuated quartz tubes, followed 
by cooling to room temperature over three weeks. The 
ordering of the chalcogen layers in our Bi 2 Te2Se single 
crystals is confirmed by the X-ray powder diffraction pat- 
terns by comparing those from Bi 2 Te2Se and Bi 2 Te3 as 
shown in Fig. 1(b), where the characteristic peaks, which 
are intensified in the ordered Bi2Te2Se compound^ are 
indicated by arrows. Note that the ordering does not 
cause a doubling nor a symmetry change of the unit cell. 

For transport measurements, cleaved crystals were 
aligned using the X-ray Laue analysis and cut along the 
principal axes, and the (111) surface was protected by de- 
positing an alumina thin film after cleaning the surface 
by bias sputtering with Ar ions. Ohmic contacts were 
prepared by using room-temperature cured silver paste. 
The sample reported here was 0.51-mm wide and 0.26- 
mm thick, with the voltage contact distance of 0.86 mm. 
The resistivity p xx and the Hall resistivity p yx were mea- 
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FIG. 1: (Color online) (a) Layered crystal structure of 
Bi2Te2Se showing the ordering of Te and Se atoms, (b) Com- 
parison of the X-ray powder diffraction patterns of Bi2Te2Se 
and Bi2Te3. Arrows indicate the peaks characteristic of 
Bi2Te2Se. (c) Temperature dependence of p xx ; inset shows 
the Arrhenius plot, (d) Plot of the conductivity a xx (= p xx ) 
vs. T -1/4 . Dashed line is the fitting of the 3D VRH behav- 
ior a xx ~ exp[— (T/To) -1 / 4 ] to the data; deviation from the 
fitting at low temperature, shown by solid line, signifies the 
parallel metallic conduction. 

sured simultaneously by a standard six-probe method by 
sweeping the magnetic field between ±14 T, during which 
the temperature was stabilized to within ±5 mK. The 
sweep rate was 0.3 T/min. 

Figure 1(c) shows the temperature dependence of the 
resistivity p xx of this ordered compound. We observe al- 
most two orders of magnitude increase in p xx upon cool- 
ing from room temperature, which is an indication of an 
insulating behavior. Indeed, the Arrhenius plot [inset of 
Fig. 1(c)] shows an activated temperature dependence 
in the range from 300 K down to ^150 K with an ex- 
citation energy A of about 23 meV. Below ^150 K, the 
transport is understood as a parallel circuit of an insu- 
lating component characterized by a 3D VRH behavior— 
and a metallic component that saturates below 10 K [Fig. 
1(d)]. 

The saturation of the resistivity at low temperature im- 
plies a finite metallic conductivity at T — K. To clarify 
the nature of this metallic state, we employed the SdH 
oscillations, whose angular dependence in tilted magnetic 
fields can provide the information about the size and the 
shape of the Fermi surface (FS) and, more importantly, 
about the dimensionality of the FS. In our Bi2Te2Se crys- 



tals, we observed SdH oscillations in both p xx and p yx 
(the latter presenting more pronounced oscillations) , and 
the two show essentially the same frequency with a phase 
shift of approximately ttJ^- Figure 2(a) shows the deriva- 
tive of p yx with respect to the magnetic field B plotted 
in 1/B± = 1/(B cos9) coordinates, where 9 is the angle 
between B and the C3 axis as shown schematically in 
the bottom inset. Two important features can be read- 
ily recognized: First, dp yx /dB is periodic in the inverse 
magnetic field, indicating the existence of a well-defined 
FS. Second, oscillations depend solely on B±, implying 
a 2D character. Note that the oscillations quickly disap- 
pear with increasing 9 (above ^40° they are hardly dis- 
tinguishable) because the amplitude of the oscillations 
strongly depends on the magnetic-field strength. The 
signature of a 2D FS can be also seen in the plot of the 
oscillation frequency vs. 9, which shows the characteris- 
tic 1/ cos 9 dependence [upper inset of Fig. 2(a)]. 

The obtained frequency for 9 ~ 0° of 64 T, which is 
directly related to the Fermi-surface cross section A via 
the Onsager relation F = (hc/2ire)A, gives &f = 4.4x 10 6 
cm , which means the surface charge-carrier concentra- 
tion N s = kp/Ai: — 1.5 xlO 12 cm~ 2 for a spin-filtered 
surface state. It is important to notice that our mea- 
sured kp is too large if the oscillations come from the 
bulk; for example, the carrier concentration for a 3D el- 
lipsoidal FS that might be consistent with our SdH data 
would be > lxlO 19 cm" 3 , which is orders of magnitude 
larger than what we obtain from the Hall data described 
later. 

Fitting the standard Lifshitz-Kosevich theory^ to the 
temperature dependence of the SdH amplitudes [inset 
of Fig. 2(b)] gives the cyclotron mass m c = 0.11m e , 
where m e is the free electron mass. Assuming that the 
electrons are Dirac-like, one obtains the Fermi velocity 
VF=hkF juic = 4.6x10 cm/s. This vp is consistent with 
the ARPES dataj 16 affirming the Dirac-fermion assump- 
tion. The position of the surface Fermi level Ep measured 
from the Dirac point can be estimated from vf and kp to 
be 130 meV, suggesting that the observed surface carriers 
are electrons. 

Once m c is known, the Dingle analysis [shown in Fig. 
2(b)] uncovers the scattering time r through the Din- 
gle temperature Td [= ^/(27rr/cs)]; the slope of the lin- 
ear fit to the data yields Tr> of 25.5 K, which gives t = 
4.8xl0~ 14 s. Hence, one obtains the mean free path on 
the surface ^f dH = vft ~ 22 nm and the surface mo- 
bility p s s dii = (e£ s s dii )/{hk F ) « 760 cm/Vs. Note that 
both quantities are underestimated, because r obtained 
from the SdH effect reflects scattering events in all direc- 
tions, whereas in the transport properties the backward 
scattering, which is prohibited in topological insulators, 3 
plays the most important role. 

In the SdH oscillations, the resistivity oscillates as 
Apxx ~ cos[2-7r(-g + i + j8)], where 27r/3 is the Berry 
phased It is knowr*i2r— that the cyclotron orbit of an 
electron in a magnetic field acquires a Berry phase n if 
its energy dispersion is linear near the degenerate point 
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FIG. 2: (Color online) (a) dp yx /dB vs. 1/B± [= l/(Bcos0)] 
in magnetic fields tilted from the C3 axis at different angles 8, 
where all curves are shifted for clarity; lower inset shows the 
schematics of the experiment and the definition of 8. Equidis- 
tant dashed lines in the main panel emphasize that the posi- 
tions of maxima for any 8 depends only on B±; upper inset 
shows the 1 / cos 8 dependence of the oscillation frequency. 
Both point to the 2D FS. (b) Dingle plot for the oscillations 
in Ap yx , which is obtained after subtracting a smooth back- 
ground from p yx , giving Td = 25.5 K; inset shows the T 
dependence of the SdH amplitude for 8 ~ 0°, yielding m c 
= 0.11m e . (c) Landau- level fan diagram for oscillations in 
dpxx/dB measured at T = 1.6 K and 8 ~ 0°. Inset shows 
dpxx/dB vs. 1/B after subtracting a smooth background. 
Minima and maxima in dp xx /dB correspond to n + j and 
ri+f, respectively. The linear fit intersects the axis at n — 
0.22. (d) Temperature dependence of the low-field Rh; dotted 
line represents the activated behavior. Inset shows the p y x{B) 
curve at 1.6 K and its fittings with the two-band model. 



(called Dirac point), whereas in ordinary metals /3 = 
0. Experimentally, the phase of the oscillations can be 
obtained from the Landau-level fan diagram, which is 
shown in Fig. 2(c) for the oscillations in dp xx /dB~2- The 
values of 1/B corresponding to minima and maxima in 
dpxx/dB [shown in Fig. 2(c) inset] are plotted as a func- 
tion of the Landau level number n<21 The extrapolation 
of the linear fit to the data shown in Fig. 2(c) suggests 
ft = 0.22; however, since our data are far from the ori- 
gin, this extrapolation involves a large uncertainty and 
we can only conclude that our data suggest a finite Berry 
phase and are not inconsistent with the Dirac nature of 
the surface state. We note that a clear it Berry phase has 
never been observed in Trs ^ 10 ' 12 

To elucidate the bulk contribution to the transport 
properties, the Hall data is useful. As shown in Fig. 
2(d), the low-field Hall coefficient Rh (= p yx /B near 



B = T) changes sign from positive to negative upon 
cooling, signifying the change of dominant charge carries 
from holes to electrons. At high temperature (> 150 K), 
the behavior of Rh (T) is thermally activated; the dotted 
line is an Arrhenius-law fitting, which signifies the acti- 
vation of holes with an effective activation energy A* = 
33 meV. This A* is of the same order as we found for 
the resistivity [inset of Fig. 1(c)], and the small differ- 
ence is attributed to the temperature dependence of the 
mobility. The prefactor of the activated behavior gives 
an estimate of the acceptor concentration N a ~ 9xl0 18 

When we look at the magnetic-field dependence of py X 
at 1.6 K, the low- field and high-field slopes are essen- 
tially different [inset of Fig. 2(d)]. The low-field Rh is 
—200 cm 3 /C, giving the effective carrier concentration of 
3.1 xlO 16 cm~ 3 . As we know from SdH oscillations, the 
surface-electron concentration is 1.5 xlO 12 cm~ 2 , corre- 
sponding to an effective 3D Hall coefficient of —5.4 x 10 4 
cm 3 /C, which is much larger than the observed value. 
Therefore, surface electrons alone cannot account for the 
low-temperature Rh and there should be other charge 
carriers in the system. Their concentration and type can 
be inferred from the high-field slope of p yx (B), since the 
high-field limit of Rh is determined only by the number 
(and type) of carriers irrespective of their mobilities; in 
our data, the high-field slope is —26 cm 3 /C, which points 
to the existence of bulk electrons with the concentration 
n of 2.4xl0 17 cm~ 3 in addition to the surface electrons. 
This example demonstrates that estimations of the bulk 
carrier density based on the low- field Rh, though often 
done in TIs, 6 ~— can be too optimistic. 

The above analysis indicates that at low temperature 
the surface electrons and bulk electrons are contribut- 
ing in parallel. It turns out that the standard two-band 
model 2 ^ 
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fits the whole p yx (B) curve very well [solid line in the 
inset of Fig. 2(d)]; here, R n and p n are the Hall coeffi- 
cient and resistivity of the bulk electrons, R s — t/(eN s ) 
and p s = pat, with the surface sheet resistance and 
t the sample thickness. This fitting yields the surface 
mobility p s = 1450 cm 2 /Vs and the bulk mobility p n 
= 11 cm 2 /Vs, along with N s = 1.5xl0 12 cm~ 2 and n 
— 2.4xl0 17 cm~ 3 . The obtained p s is about two times 
larger than the /itf dH estimated from the SdH analysis 
[Fig. 2(b)]. This is expected for the topological sur- 
face state as we mentioned before. On the other hand, 
the bulk mobility p n of 11 cm 2 /Vs is near the boundary 
between the band and hopping transport regimes ^ sug- 
gesting that the bulk electrons move in a very disordered 
potential. The fraction of the surface contribution in the 
total conductance at 1.6 K is calculated to be ~6%, which 
is about 20 times larger than that achieved^ in Bi2Te3. 

From the detailed information we have gathered, a 
comprehensive picture emerges for the bulk transport 
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FIG. 3: (Color online) Schematic picture of the bulk and sur- 
face band structures (left), together with the energy diagram 
of the density of states of the bulk and impurity bands (right). 
In the impurity band which is due to the acceptor levels and 
is located within the band gap, only the central part forms 
the extended states and the tails consist of localized states. 

mechanism in Bi2Te2Se: Taking into account the rela- 
tively large concentration (9xl0 18 cm~ 3 ) of acceptors, 
it is reasonable to assume that those acceptors form an 
impurity band (IB) within the energy gap of Bi2Te2Se; 
note that the IB is formed because the wave functions of 
the electrons bound to impurity sites overlap to form ex- 
tended states, which are responsible for the finite n-type 
bulk conduction at zero temperature with n — 2.4xl0 17 
cm~ 3 . The chemical potential is pinned to this IB at low 
temperature and is obviously located within the bulk en- 



ergy gap (Fig. 3). As temperature increases, hopping 
conduction of electrons using localized states becomes 
possible, giving rise to the VRH behavior which occurs 
in parallel with the degenerate IB conduction. At even 
higher temperature, thermal activation of electrons from 
the bulk valence band to the acceptor levels takes place, 
leading to the p-type bulk conduction. 

In conclusion, we have shown that Bi2Te2Se is the best 
material to date for studying the surface quantum trans- 
port in a TI. The surface contribution in the total con- 
ductance of our Bi2Te2Se crystal is ~6%, which is the 
largest ever achieved in a TI. The bulk mobility of 11 
cm 2 /Vs indicates that the IB conduction in our sample 
is at the verge of localization, and not much further re- 
duction in the number of acceptors would be needed to 
quench the degenerate bulk conduction. Once such crys- 
tals become available, they will allow us to study the 
plethora of topological quantum phenomena that have 
been predicted 3 - - — for this class of materials. 
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